Mott transition in the Hubbard model on the hyper-kagome lattice 
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Motivated by recent experiment on the Na4lr308 compound we study the Hubbard model on the 
"hyper-kagome lattice", which forms a three-dimensional network of corner sharing triangles, using 
dynamical cluster approximation (DCA) method with Nc=12 combined with the continuous-time 
quantum Monte Carlo (CT QMC) method. The system undergoes a Mott transition if the Hubbard 
interaction U /W {W is the bandwidth) exceeds the value of 1.2 for T=0.1 and displays reentrant 
behavior due to competition between the magnetic correlation and the kinetic energy of electrons 
due to the geometrical frustration. We observe a "critical slowing down" of the double occupancy 
which shows evidence of a continuous transition. The nearest-neighbor and next nearest-neighbor 
spin-spin correlations indicate a paramagnetic metalhc state in the weak-coupling regime and an 
antiferromagnetic (AF) Mott insulator in the strong-coupling regime within the temperature range 
which we can access with our numerical tools. 

PACS numbers: 71.10.Fd 



Frustrated electronic systems exhibit a wide range of 
phases due to competition between strong electronic cor- 
relations and geometric frustration. The Na4lr308 com- 
pound with a S=l/2 three-dimensional network of cor- 
ner sharing triangles, which is called the "hyper-kagome 
lattice", was discoveredji Due to the structure of cor- 
ner sharing triangles, the system displays the electronic 
correlations and geometric frustration at the same time. 
Each tetrahedron in the pyrochlore lattice, which is 
shown in Fig. 1, consists of Ir (indicated as "filled" 
circles) and Na (indicated as "empty" circles). S=l/2 
spins are carried through Ir^+ on the hyper-kagome lat- 
tice (blue lines in Fig. 1) because Ir is tetravalent with 
five electrons in 5(i orbitals. In recent experiment by Y. 
Okamoto (Ref. 1), a spin liquid state at low temperature 
was found. Theoretical work, related to this experiment, 
has mostly addressed the strong-coupling limit using the 
Heisenberg model^'^'^'^'^. These papers concentrated on 
the nature of the ground state which was identified as 
the spin liquid state. Even more recently, it was shown 
experimentally that this material undergoes a transition 
to a metallic state under pressure^. 

The two important theoretical questions which we will 
address in this Letter are if the phase diagram shows a 
reentrant Mott phase with a first-order transition as was 
found for the anisotropic triangular lattice-!^ and if the 
spin liquid state is stable at finite temperatures. The 
dynamical mean field theory (DMFT)^!^ is a success- 
ful method to describe the Mott transition. However, 
since the single-site DMFT ignores the spatial fluctua- 
tions, the system becomes metallic as temperature is de- 
creased. This behavior is quite different from the phase 
diagram obtained from quantum Monte Carlo (QMC) on 
the triangular lattice. ^"'^^ The cluster extension of the 
DMFT methods such as the dynamical cluster approx- 
imation (DCA)i^ii2iii^ and the cellular dynamical mean 
field theory (CDMFT)^^ were developed in order to con- 
sider the spatial fluctuations. In these methods the short 
range correlations are treated exactly within the cluster 
size, while the long range correlations are included on the 




FIG. 1: In the pyrochlore lattice, "filled" circles and "empty" 
circles present Ir and Na, respectively. The blue fines show the 
hyper-kagome lattice of Ir. a is the lattice constant. Lattice 
points labeled by 1,2,.. .,12 represent relative positions in each 
unit cell. 



mean field level only. The difference between both meth- 
ods lies in the boundary condition of the cluster given 
by the Laue function^'^. In this paper, we employ the 
DCA method combined with the continuous-time quan- 
tum Monte Carlo (CT QMC) methodic as cluster solver 
and measure the total density of states (DOS), double 
occupancy, nearest-neighbor and next nearest-neighbor 
spin-spin correlation functions. 

We consider the Hubbard model on the hyper-kagome 
lattice at half-filling: 



(1) 



where Cia-(cJ^) is the annihilation (creation) operator of 
an electron with spin a at the i-ih site, t is the hopping 
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matrix element and U represents the Coulomb repulsion. 
We set t=1.0 and the bandwidth is W^6t. For the DCA 
method we consider an unit cell which has twelve sites 
(A'c=12) in the simple cubic lattice, as shown in Fig. 1. 
In the DCA method the Green's functions are determined 
by the self-consistency condition. We directly write the 
hopping matrix T(k) with the periodic boundary condi- 
tion in the real space. Here Gq, S and T(k) are described 
by 12 X 12 matrices. In this case the self-consistent equa- 
tion for the DCA method is the same as for the CDMFT 
method except the hopping matrix T(k): 



dk- 



S(»Wn), (2) 



where /x is the chemical potential and w„ is the Matsub- 
ara frequency. The Brillouin zone of the superlattice also 
becomes simple cubic lattice with < k^, ky, kz < 7r/a, 
where a = 1 is the lattice constant and the summation 
of k is taken over the Brillouin zone of the superlattice. 
We use 8 x 10^ k-points for integration of Eq. We 
solve the quantum impurity problem using the CT QMC 
method. Unlike the Hirsch-Fye discrete time quantum 
Monte Carlo methodic which is carried out by the aux- 
iliary Ising-like spins, the CT QMC method performs a 
random walk in the space of the perturbation expansion 
without Trotter decomposition error. The computational 
time of CT QMC scales as (k) ~ O.bNcUfJ and is there- 
fore superior to the Hirsch-Fye quantum Monte Carlo 
method for which the computational time scales as (k) ~ 
2NcU(3, where (fc) is the matrix size. We briefly summa- 
rize the CT QMC method. The main idea of CT QMC is 
to perform a diagram expansion of the partition function 
Z = 2'oe~'^-^'''^"T"i in powers of the interaction U. In 
this case we can reexpress the partition function as 

^ kl J i ^ (3) 

< n|(Ti)n|(Ti) • ■ •n|(Tfc)n|(Tfe) >, 

where < n|(Ti)n|(Ti) • ■ • n-f (rfc)ri|(rfe) > is determined 
by the non-interacting Green's function and Wick's the- 
orem and Zf) is the unperturbed part. The interacting 
Green's functions are calculated by numerical averaging 
of Eq. (3). We use 3 x 10^ QMC sweeps approximately 
and consider twenty-one Matsubara frequency points. 

Now let us investigate the phase diagram of the hyper- 
kagome lattice Hubbard model at half filling. In order to 
check the existence of a gap, we calculate the density of 
state (DOS) at the Fermi level (w = 0) by using the Fade 
approximation and the analytic form which is given by 

G(/3/2) « ^p{lo = 0), (4) 

where P is the inverse temperature and p(ti^) is the den- 
sity of states. Because Eq. (4) is satisfied at suffi- 
ciently low temperature, we employ the Fade approxi- 
mation in T>0.1 regions and compare the results of both 
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FIG. 2: Phase diagram of the Hubbard model on the hyper- 
kagome lattice 
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FIG. 3: The density of state corresponding to (a) C//W=1.17 
and (b) [//VF=1.25 for T=0.1. The Fade approximation is 
employed for analytical continuation 



methods at T=0.1. The Fade approximation method 
for analytical continuation works well at the high tem- 
perature regime due to very small error in CT QMC 
method. Moreover, the results of Eq. (4) support those 
obtained from the Fade approximation at low tempera- 
ture. The phase diagram in U-T plane is shown in Fig. 
2. As temperature is decreased in the transition region 
(f7/VF=1.28) above T=0.3, the system prefers the metal- 
lic state due to the dominant low energy excitation of 
electrons. Below T=0.2 the system is strongly controlled 
not by itinerant electrons but by antiferromagnetic (AF) 
fluctuations, which are ignored in the single-site DMFT. 
Therefore, the Mott insulator lies in the low temperature 
regions. This results in a reentrant behavior which is 
also presented in the Hubbard model on the anisotropic 
triangular lattice^^ related to the frustrated organic ma- 
terial K-(BEDT-TTF)2Cu[N(CN)2]CliSiS0. We analyze 
the DOS close to critical U jW at T=0.1 using Fade ap- 
proximation for analytical continuation. The DOS for 
UIW=\.\1 and UlW=\.2h is exhibited in Figs. 3(a)- 
(b). At U/W=1.17 a quasiparticle peak due to itinerant 
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FIG. 4: Double occupancy as a function T for (a) U /W=1.Q, 
(b) U/W=l.2 and (c) U/W=1.2>2,2>. 



of electrons is formed strongly around the Fermi level. At 
U/W=l.25 the DOS changes dramatically at the Fermi 
level and the system turns to a Mott insulator. We did 
not observe any evidence for the pseudogap formation 
like the results of kagome lattice^. The magnetic state 
is suppressed due to the geometrical frustration in con- 
trast to the square lattice^^i. To investigate the character 
of Mott transition in more detail we present the double 
occupancy as a function of temperature for various in- 
teractions in Figs. 4(a)-(c). As temperature is decreased 
at weak interaction strength U/W=1.0 in Fig. 4(a), the 
double occupancy is increased. This means that the en- 
tropy due to kinetic energy of the dominant quasiparti- 
cles, is optimized by pushing the interaction in order to 
decrease the free energy. Non-monotonic dependence of 
the double occupancy is observed at intermediate inter- 
action U/W=1.2 in Fig. 4(b). For T >0.3 the system 
displays the insulating state due to effect of the dominant 
local magnetic moment. As temperature is decreased, 
it changes into the metallic state due to development 
of the quasiparticle. After passing the peak tempera- 
ture To=0.15, at lower temperatures the system reenters 
the insulating state because the magnetic correlations 
get enhanced. This is indeed evidence for the reentrant 
behavior which is exhibited in geometrically frustrated 
systems like the case of Hubbard model on anisotropic 
triangular lattice. In Fig. 4(c), we see a small gradual 
increase in the double occupancy as temperature is de- 
creased, reflecting a gain in kinetic energy. Our analysis 
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FIG. 5: Double occupancy as a function U /W for T=0.2 
and r=0.1. There is a "critical slowing down" behavior at 
r=0.1 which represents the evidence of a continuous transi- 
tion. Uc/W^l.2 for r=Q.l 



in the Fermi liquid region {U/W— 1.0) and Mott insula- 
tor region {U/W=1.33) is similar to that obtained from 
the cluster-extension of DMFT method on the square 
lattic^. We further investigate the double occupancy 
as function of interaction at various temperatures in or- 
der to see the type of transition. The result is shown in 
Fig. 5. As interaction is increased, at T—Q.2 the double 
occupancy decreases smoothly because of the develop- 
ment of local magnetic moment. At T=0.1 the "critical 
slowing down" behavior is observed around Uc/W=1.2. 
This behavior is quite similar to the results of the fully 
frustrated Hubbard model which indicates a second-order 
transition in the single-site DMFT calculation^i. We sus- 
pect that there is a continuous transition in the Hub- 
bard model on the hyper-kagome lattice. We calculate 
the nearest-neighbor {SfSfj^^) and next nearest-neighbor 
spin-spin correlation (SfSf^2)j which are shown in Fig 6. 
In the weak interaction regions the system displays para- 
magnetic metallic behavior, unlike the case of the square 
lattice which has strong AF correlations due to perfect 
nesting. As interaction is increased, there is a compe- 
tition between the quasiparticle formation and the frus- 
trated spin correlation. At Uc/W=1.27 for T=0.2 both 
the nearest-neighbor and next nearest-neighbor spin-spin 
correlations are enhanced at the same time. We find the 
AF Mott insulator in the strong interaction regions. We 
do not find a spin liquid state for r>0.1. 

Finally, we compare our results with those recently 
obtained from the mean-field calculation^^. In the 
mean-field calculation the system has continuous metal- 
insulator transition at critical interaction Uc/W=1.03. 
These results are comparable to our results with a contin- 
uous metal-insulator transition at Uc/W=1.2 for T=0.1. 

In summary, we have presented results of DCA study 
combined with CT QMC method as a cluster solver in 
the Hubbard model on the hyper-kagome lattice. The 
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FIG. 6: The nearest-neighbor and next nearest-neighbor 
spin-spin correlation function as a function of U/W for 
r=0.2. Both spin-spin correlations are rapidly increased at 
Uc/W=1.27 



phase diagram in the U-T plane has reentrant behav- 
ior due to geometrical frustration and is quite similar 
to that on the anisotropic triangular lattice. The differ- 
ence of both lattice types is that while the hyper-kagome 
lattice exhibits a continuous transition, the anisotropic 
triangular lattice displays a first-order transition. Us- 



ing Padc approximation for the analytical continuation 
we calculated the DOS close to critical Uc/W=1.2 for 
r=0.1. There is rapid change in the DOS at Fermi level. 
The formation of pseudogap is not observed because the 
magnetic fluctuations are suppressed. We investigated 
the double occupancy as a function of temperature and 
as a function of interaction U/W to study the nature 
of the Mott transition in more detail. For intermediate 
values of the interaction, the double occupancy as a func- 
tion of temperature has non-monotonic behavior which 
represents the character of reentrant behavior. We also 
find the "critical slowing down" of the double occupancy 
as a function of interaction. We suspect that at low tem- 
perature there might be a signal of a continuous tran- 
sition. Finally, we calculated the nearest-neighbor and 
next nearest-neighbor spin-spin correlations as a func- 
tion of interaction. We observed a paramagnetic metal- 
lic state in the weak-coupling regime. When the interac- 
tion passes the critical point, both nearest-neighbor and 
next nearest-neighbor spin-spin correlations are rapidly 
increasing. In the strong-coupling regions the system 
shows the AF Mott insulator state. We did not observe a 
spin liquid state (paramagnetic Mott insulator) at lower 
temperature. However, from our calculation we expect 
that a continuous Mott transition with reentrant behav- 
ior will be observed experimentally. 

We would like to thank A. J. Millis for helpful discus- 
sions. 
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